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We present a unified study of the effect of periodic, quasiperiodic and disordered potentials on 
topological phases that are characterized by Majorana end modes in ID p-wave superconducting 
systems. We define a topological invariant derived from the equations of motion for Majorana 
modes and, as our first application, employ it to characterize the phase diagram for simple periodic 
structures. Our general result is a relation between the topological invariant and the normal state 
localization length. This link allows us to leverage the considerable literature on localization physics 
and obtain the topological phase diagrams and their salient features for quasiperiodic and disordered 
systems for the entire region of parameter space. 
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Introduction - Recent claims of the detection of Ma- 
jorana fermions in semiconducting/superconducting het- 
erostructures have stirred new excitement leading to sev- 
eral avenues of inquiry [TJ [2] . A major concern in these 
effectively spinless p-wave superconducting wires is the 
role of spatially varying potentials, be they externally 
applied or due to disorder. The latter has in fact been 
a rich, active topic of study for over a decade in terms 
of delocalization-localization physics in one-dimensional 
(ID) class D systems [21 E]. The pioneering work of 
Ref. [3] specifically explored the conditions for the ex- 
istence of Majorana boundary modes in these disordered 
systems, arguing that a finite amount of superconductiv- 
ity is required to drive the system into such a topological 
phase. More recently, several works have further inves- 
tigated this aspect, primarily exploring either weak (or 
slowly varying) disorder or a special point in parame- 
ter space corresponding to the disordered quantum Ising 
chain [5j E] - In this Letter, we perform a cohesive study 
of the topological phase diagram for a range of poten- 
tial landscapes on a lattice, extending and unifying work 
on periodic [7J [5], incommensurate [5], and disordered 
potentials HU| fTT]. 

Our central observation links the normal state localiza- 
tion properties to the strength of the superconducting 
pairing required to engender an end Majorana mode, thus 
enabling us to obtain the topological phase diagram for a 
variety of situations (see Fig. [TJ . We define a topological 
phase (T-phase) to be one in which there are zero-energy 
Majorana modes at the ends of an infinitely long sys- 
tem. Specifically, we define a topological invariant based 
on the decay length of the end mode [7] and provide 
a relationship between this length and the correspond- 
ing zero-energy localization length of the normal system. 
Unlike in the uniform case, we find that even in the sim- 
plest case of periodic potentials, the T-phase requires a 
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FIG. 1: Topological phase diagrams (T - topological; N - 
non-topological) as a function of superconducting gap, A, 
and potential strength, W, for (a) a periodic potential having 
the pattern (W, W, — W, — W), (b) a quasiperiodic potential 
W cos(27rum) for any irrational u), (c) box-distributed (uni- 
form) disorder, and (d) Lorentzian distributed disorder. 

finite amount of superconductivity, which thus acts as 
a new knob to access the phase. In the quasiperiodic 
case, the topological boundaries reflect the morphology 
of the fractal patterns exhibited by the normal state en- 
ergy spectrum (see Fig.[2|. 

For disordered systems, our observation proves pow- 
erful in that it allows us to leverage the vast body of 
literature on normal state localization physics to identify 
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topological properties of the disordered superconductor. 
The topological regions arise as a result of the competi- 
tion between the localizing effects of disorder and the su- 
perconducting pairing which tends to spatially separate 
the Majorana fermions composing a Dirac state. Our 
analyses provide a mapping for the phase boundary be- 
tween strong and weak disorder limits and reproduce the 
exact form for the intermediate point corresponding to 
the random field transverse Ising model. We present two 
representative cases of our analyses for the disordered 
phase diagram in Fig. [ljc-d) . 

Model and Topological features - Turning to our start- 
ing point, we model the superconducting wire as a ID 
tight-binding system of spinless electrons exhibiting p- 
wave superconductivity described by the TV-site Hamil- 
tonian 12J 

M 

H = E [ - * (tif n+l + 4+l/n) + *(/»/„+! 
n=l 

+ 4+i/J)-M4/»-i/2)], (!) 

where t is the nearest-neighbor hopping amplitude, A 
is the superconducting gap function (taken to be real), 
and we will eventually take TV — > oo. The various 
cases of periodic, quasiperiodic, and disordered poten- 
tials are encoded in the local on-site chemical potential 
\x n and are characterized by a typical potential strength 
W. The Dirac fermion /„ can be expressed in terms of 
Majorana fermions, f n = (a n + ib n )/2, which are Her- 
mitian operators satisfying the anticommutation rules 
{a n ,a m } = {b n ,b m } = 2S n<m and {a n ,b m } = 0. 

We now construct a topological invariant that links 
topology to the eigenvalue structure of zero-energy 
boundary modes described by the Hamiltonian in Eq. (JlJ . 
Specifically, the end Majorana modes that decay into 
the bulk can be represented by the operators Q a — 
J2n a n a n, Qb = J2 n ^nb n , where the wave function a n 
obey the zero-energy equations of motion derived from 
Eq. 0. As detailed in Ref. [7J, these equations can be 
represented in the transfer matrix form 

¥ )■ 

(2) 

Since the A n may be taken as functions of fJ, n /t an d A/i, 
we set t = 1. A similar expression holds for the /3 n . 

The existence of end Majorana modes requires the a n 
(or P n ) to be normalizable. We denote the number of 
eigenvalues of the transfer matrix A(W, A) = n^Li A n 
with magnitude less than 1 by rif. For rif = 0,2, the 
system is topological, whereas for rif = 1 it is non- 
topological [7 . Since the topology of the system depends 
only on the magnitude of A, we take A to be positive; 
hence | det A\ < 1. Then the two eigenvalues of A obey 
| Ai A 2 1 < 1. Therefore, for |Ai| < |A 2 |, we have |Ai| < 1 



TABLE I: Criteria for topological phase for a selection of pe- 
riodic potentials. 
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pattern of \x n 


topological for 
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...,W,W,W,... 


\W\ < 2 


2 


...,W, -W, . . . 


A > \W\/2 


4 


...,w,w,w,-w,... 


A 2 > W 2 /2 - 1 


4 


...,w,w,-w,-w,... 


A > W 2 /4 



and rif is completely determined by the larger eigenvalue 
A 2 . We thus define a topological invariant 

v = - (-1)"' = sgn(ln|A 2 |), (3) 

which can also be expressed as v — — sgn (/(l)/(— 1 )), 
where f(z) = det (I — Az) is the characteristic polyno- 
mial of A [7J. From the conditions on rif, we see that 
v = — 1 for the T-phase and v = 1 for the non-topological 
phase (N-phasc). 

Uniform and Periodic Potentials - As a simple appli- 
cation of our topological invariant, we consider periodic 
patterns in the sign of the chemical potential (detailed in 
Ref. [7J). Table I presents the conditions for T-phases for 
some select patterns. The comparison of the phase dia- 
grams for the uniform case and simple periodic potentials 
highlights unusual aspects of the former's phase diagram. 
In particular, the system is topological for \W/t\ < 2 for 
any A / 0. This stems from the fact that the system 
has a bulk gap at E = as long as W lies in the range 
[— 2t, 2t] [T5]. In contrast, non-uniform potentials tend 
to open a bulk gap at E = with a size that grows with 
increasing W. This leads to a phase boundary which 
generically has > 13 . We now quantify this ob- 
servation and extend it to other potential landscapes. 

Features of the Topological Phase Diagram - As ob- 
served in Ref. [1] , the product of transfer matrices which 
appears in A is strongly reminiscent of that used to de- 
termine localization properties of the normal state An- 
derson disorder problem. We build on this observation by 
determining the phase boundary from the normal state 
properties of the system. In general, this leads to a crit- 
ical amount of superconductivity required to drive the 
system into a T-phasc. 

To this end, for < A < 1, we perform a sim- 
ilarity transformation A n — y/SSAnS^ 1 with S — 
diag(<5 1 / 4 ,l/(5 1 / 4 ) and S = j^f. The matrices A n are 
of the form shown in Eq. pi) with A — > and fj, n — > 
lx n /\/l — A 2 . This immediately gives 

(4) 

Taking the logarithm of the eigenvalues of Eq. Q, the 



condition that | A2 1 = 1 is given by 
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1. A + A 
= 2 ln T^A 



(5) 



where we have defined the Lyapunov exponent of the 
normal state system, 7(W) = lim^^ -h In |Aa(Wi 0)|. 
Eq. ([5]) thus describes the phase boundary separating the 
topologically trivial and non-trivial regions of the phase 
diagram. This relation quantifies the observation in [4j 
110] that in general a critical amount of superconductivity 
must be applied before the system is driven into a T- 
phase. For the case in which the system is metallic (i.e., 
j(W) = 0), any non-zero A will give rise to a T-phase. 

The form of the phase diagram for A > 1 may be 
obtained by noting that the transformation 

M„ Mn/A, A ->• 1/A , and P -> P (6) 
where P ^ P : { M „} -> {(-l)Vn}, 

leaves the eigenvalues of A unchanged for AT even. Thus, 
if a point (Wq,Ao < 1) lies on the phase boundary of 
P, then (Wo/Ao, 1/A ) lies on the phase boundary of P. 
This duality strongly constrains the form of the phase 
boundary in the cases where the distribution is invariant 
under the transformation in Eq. Q. 

Finally, at the point A = 1, the system maps to the 
well-studied quantum Ising chain subject to a random 
transverse field [5J [6] . The matrix ^(W 7 , 1) has the eigen- 
values rinLi an d 0. Eq. reveals that the phase 
boundary passes through the point for which 



(In |^„|) =ln2, 



(7) 



where (ln|jU„|) = lini/v^oo -fo J2 
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1 ]n\fj, n \. These rela- 
tions allow us to obtain the superconducting topological 
phase diagram for quasiperiodic and disordered poten- 
tials purely based on the normal state localization prop- 
erties. 

Quasiperiodic Potentials - Here we consider two cases 
of potentials to study periodicity that is incommensurate 
with the underlying lattice. In the first instance, where 
M« = W cos (27rwn) and uj is irrational, the normal state 
features a well-studied metal-insulator transition at the 
critical value of W = 2 [14] . The normal state Lyapunov 
exponent takes the form Jqp(W) = In (W/2) for W > 2 
and for < W < 2 for u irrational QH \T§\. Eq. ^ 
then predicts a T-phase for 




FIG. 2: (a) The Lyapunov exponent yn of the normal state 
(A = 0) for /i n = V + 2 cos (2-Kun) with uj = 1/10. (b) Topo- 
logical phase boundary showing the merging of the topological 
regions (T) as described in the text. For A < 1, there are are 
10 distinct regions of non-trivial topology (T). At A = 0.2, 
the four central regions have merged to form a single region, 
(c) A colorscale plot of yn as a function of V and uj. Darker 
regions correspond to lower values of jh- Striations show 
the spectrum's sensitivity to q for uj — p/q. (d) The topo- 
logical phase diagram (A = 0.2) reflects the geometry of the 
low- lying values of yn in (c). 



A second example of particular interest is the Harper 
potential \i n = V + 2 cos {2-kujti), corresponding to the 
problem of an electron hopping on a 2D square lattice 
with each plaquette enclosing a magnetic flux [17 . The 
associated normal state energy spectrum, the celebrated 
Hofstadter's butterfly, has a rich fractal structure stem- 
ming from the fact that w can take irrational values and 
that values of oj = p/q (p and q relatively prime) give rise 
to q bands separated by q — 1 non-zero gaps (Fig. J2ja)). 
From Eq. Q , we directly obtain the topological phase di- 
agram for non-zero A. As shown in Fig.[2]jb), for A <C 1, 
there are q topological regions inherited from the nor- 
mal state which fuse as A is increased. In the U) — V 
parameter space, as expected from our general analysis, 
the normal state properties (Fig. ffic)) directly inform 
the topological phase diagram (Fig. |2[d)). 



A>\ W - 



1. 



(8) 



This result holds for all values of A given that the trans- 
formation lo — > uj + 1/2 yields Eq. ([7| and that the du- 
ality transformation, A —> 1/A and W — > W/ A, leaves 
Eq. ([8| invariant. Finally, Eq. ([7| also shows that the 
point (W, A) = (4, 1) lies on the phase boundary. 



Disordered Potentials - We begin with our most gen- 
eral results for the disordered topological phase diagram 
which pertain to the limits of weak and strong disor- 
der. Consider weak, uncorrelated disorder satisfying 
(MnM«'> = U5 njn >, = 0. Using the known Lya- 

punov exponent obtained from perturbation theory for 
the normal state system [T5J \T§\ and Eq. ^ , we obtain 
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the condition for the T-phase 



A > 



f r(3/4) Y 
\T(1/4)J 



U w 0.114 C/. 



(9) 



This result may be compared to that of a continuum 
model based on the Dirac equation, which gives a topo- 
logically non-trivial phase for A > |f/ = 0.125J7 
(see [TO])- For disorder distributions that are symmet- 
ric around 0, the self-duality condition P = P in Eq. ^ 
is satisfied. In this case, we can employ this duality trans- 
formation, A — > 1/A and U — > U/A 2 , to show that 
Eq. ([9]) also describes the phase boundary in the limit 
of strong disorder. 

As the most generic representative for disorder, we 
now turn to the case of 'box' disorder (B) for which 
the probability of /i„ falling at any point in the range 
— W/2 < /i„ < W/2 is equally likely. The low-energy be- 
havior as shown in the numerical simulation in Fig. [l|a) 
is in good agreement with Eq. ([9| (for box disorder, 
U = W 2 /12). Eq. Q reveals that the phase bound- 
ary passes through the point (W, A) = (W c ,l), where 
W c = 4e m 10.873 (box disorder) with e being the base 
of the natural logarithm. 

A noteworthy find is the observed discontinuity suf- 
fered by the phase boundary as it passes through the 
random field quantum Ising point A = 1 (Fig. fljc-d)). 

To understand its origin, we note that 7 (yj=^fji the 
effective Lyapunov exponent that we seek in Eq. ([5]) cor- 
responds to that of very strong disorder for A — > 1. In 
this limit, we can use the known form of the normal state 
Lyapunov exponent for W ^> 1 [19 in Eq. (j5l), invoke 
self-duality and obtain the phase boundary to linear or- 
der around (W, A) = (4e, 1), 
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for A < 1, 
for A > 1. 



(10) 



As seen in Fig.jljc), this result is in reasonable agreement 
with numerical simulations. To go further, treating the 
quantity S = perturbatively reveals corrections to 
Eq. Q yielding (ln|/i„|) = In 2- (l + (1/ M ) 2 ) § + 0(6 2 ). 
This shows that the phase boundary is fragile towards 
singularities when [i n is allowed to come arbitrarily close 
to zero. Indeed, our simulations have shown that the dis- 
continuity is absent for disorder distributions that avoid 
zero energy. From our extensive simulations and gen- 
eral insights in the box disorder case, we conclude that a 
large class of disorder distributions that cover zero energy 
give rise to a discontinuity in the slope of the topological 
phase boundary at A = 1. 

Finally, we turn to the specific case of disorder drawn 
from a Lorentzian distribution 



fa(x;W) 



W 



1 



7T X 2 + W 2 ' 



The phase diagram is exactly soluble in this case 
since the normal state density of states is known ex- 
actly |20j . The zero-energy Lyapunov exponent, first 
obtained by Thouless [ST], takes the form 7l(IF) = 

in (w/2 + yr + W 2 /4 S j. Once again invoking Eq. (pi 

and self-duality of the phase diagram yields a phase 
boundary 



W = 



2A for A < 1, 
2 for A > 1. 



(12) 



This result, as shown in Fig. [TJc) , is in excellent agree- 
ment with numerical simulations. It should be pointed 
out that the features of this phase diagram are extremely 
unusual. For instance, Eq. ^ fails to hold because the 
second moment (i^nUn'} is ill-defined for fa. This exam- 
ple is noteworthy since, for W > 2 the system is always 
in an N-phase; no amount of A can drive the system 
into a T-phase. Studying these examples has shown us, 
among other features, that typically the larger the dis- 
order, the more superconductivity is required for Majo- 
rana end modes to exist, and that the topological phase 
diagram is highly sensitive to the nature of the disorder 
distribution. 

In conclusion, forging a connection between the normal 
state localization properties and the behavior of Majo- 
rana end modes has provided us a powerful means for 
constructing the complete topological phase diagram for 
superconducting wires. Future work would include more 
extensive utilization of known Anderson localization re- 
sults and other physical considerations, such as interac- 
tions and finite temperature. 

Note - In the very final stages of preparation of this 
manuscript, we have been made aware of another work 
which has some overlap with ours: Ref. 22J. 
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